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1. Introduction 

In this paper we study a componentwise timeshift transformation (CTT) for 
coupled dynamical systems with transmission delays, which allows one to change 
the involved delays under some restrictions. We assume a network of N dynamical 
systems whose topology is a multigraph, which is set of nodes, TV = {1, N}, and 
a multiset £ of directed edges I £ M X N '. The dynamics of the systems under 
consideration are assumed to be of the form 

(1-1) ^(^/. ((^(i-rtl)))^.) el, j = l,...,N. 

Here the functions 

s : £ — > V and t : £ — > V, 

identify source s (£) and target t (£) of a link £ e £ , and Ij — {£ £ £ : t (£) = j} is 
the input set of the node j. The time delay of the signal traveling along the link £ 
is denoted by r (£). Such a function, r : £ — > [0, oo), is called a delay distribution. 
The functions fj : M# Ij — > R describe the dynamics of the individual systems. We 
will write / = (/i, /at) T to denote the full right hand side of (jl.lj) . In general, 
boldface symbols denote vector valued variables. 

The basic idea of the CTT is to introduce new variables 

(1.2) ViW-Xjit + rij), 

by shifting the cells independently in time to some extend rjj > 0. One readily 
checks that these new variables fulfill the dynamics 

(1-3) ±Vj(t) = fj ((y sW (i-f(f))) <£j J , j = 1,...,N, 

with modified delays f (£) = r (£) — rjt^ + r] s m. 

In section|3]we explain why (jl.3l) may possess properties which make it preferable 
to (jl.ip and render the idea (|1.2|) useful. But before doing so let us pose and answer 
a question. What does knowledge about features of (|1.3[) tell us about the features 
of (ll.l[) ? It might seem intuitive that they have a lot in common since the timeshift 
Vj W ^ Uj (t ~~ Vj) inverses (|1.2p . This is true for the particular solution p.2p . but 
in general, the expression yj (t — rjj) (with rjj > 0) may not be defined at all. This 
is because, in general, solutions of delay differential equations (DDEs) cannot be 

l 
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continued backwards. Furthermore, one should introduce state spaces and flows to 
formulate and compare the dynamical properties of f)l . 1 [) and f| 1.31) . This will be 
done in a quite general setting using semidynamical systems in section [3] Before 
that we study the special cases of equilibria and periodic orbits of DDEs in section 
II 

2. Spectrum of equilibria and periodic orbits in (11.11) and (ll.3|) 

An equilibrium point x = (xi, xn) 6 Ylj=i ^™ fc °f (UTTJ) , and equally for (|1.3p . 
is a point which satisfies 

(2-1) fi((*.w) teIj ) = ' 3 = 1,-,N- 

Let us assume that fj € C 1 . Then, if their real parts are different from zero, the 
characteristic exponents of x determine its stability. A characteristic exponent of 
x in (jl.ip corresponds to an exponential solution £ (t) = e A *£ of the variational 
equation 

(2.2) ^ (i) = £ difi ((^') Wj e* (* - t (*)) . 

Obviously, the timeshifted variation x (0> given by Xj (*) = £j (^ + %) is then a 
solution of the characteristic equation of x in (11.31) 

(2-3) Xj (t) = J2 D fi (W))^ e iJ Xfc (* " •? W) • 

Hence, the characteristic exponents of x are the same in (|1.1[) and (|1.3p . 

Similarly, the stability of a periodic solution x(t) = x (t + T) is determined 
by its Floquet exponents (again, if their real parts are different from zero). Each 
exponent A corresponds to a solution £ (t) of the variational equation 

(2.4) 4- (t) = d e fj ((x fe (t t (O) W,) a- (t - t (*)) 

which has the form £ (t) = e xt p(t) with a periodic function p(t) = p(i + T). 
Again, the timeshifted solution x (t) with \j (t) = £j (t + T]j) fulfills the variational 
equation 

(2.5) X j (*) = E (* ~ f (O) Wj Xk(t~f (I)) 

leij 

of the corresponding periodic solution %jj (t) — Xj (t + r)j) in the transformed system. 
Let us summarize this section in a theorem: 

Theorem 1. (i) A fixed point x possesses the same characteristic exponents in 

rn\) and rop . 

(ii) A periodic solution x (t) possesses the same Floquet exponents in il.l]) as 
the corresponding transformed solution y (t) in il.3\) . 
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3. The Componentwise Timeshift Transformation 

3.1. The CT-transformation. In this section we set out to define the transfor- 
mation (|1.2[) in a more rigorous way. Recall that the state spaces for DDEs are 
infinite dimensional. They have to contain segments of functions, which represent 
the history of the solution x (t) . For the coupled system (|1.1[) , we choose the state 
space to be 

N 

c = n ° ([~ r ^' °i ; R " j ) ' r i = maxr w ' 

where Oj = {£ G £ : s (£) — j} is the set of all outgoing links from the j-th cell. 
Similarly, we choose 

N 
3=1 

as state space for ()1.3|) . Assuming that solutions exist for all future times (e.g. if 
fj are Lipschitz continuous), there exist semi flows 

$ : [0,oo) x C -> C, (i, a;) i->- $ t (as) , 

* : [0,oo) xC^C, (t,y)H- ¥*(»), 

for (fTTTjl and (fl~3]) . Now let us formulate (fl~2]) as a state space transformation 
T : C — > C. In a natural way we define T componentwise, for j = 1, TV, as 

(3.1) r i [x ](t) = /^ (a ' o)) i (0) > ^[-minfef^O], 

Recall that we need rjj > 0. Further we may assume min^ rjj = since f (£) is 
invariant under a simultaneous shift rjj i-> r/^ + a of all cells by a common amount 
a. Note that necessarily r^- — fj > r,-, so (|3.1[) is well-defined [see Figure [3~T] . One 
easily checks that 

(3.2) ro$ = $oT. 

That is, T transforms solutions of (|1 . 1|) into solutions of (|1.3[) . As inherited from 
the semi flow <&, the transformation T is neither injective nor surjective in general. 
Therefore, one cannot expect a dynamical equivalence of (|1.1[) and (|1.3p in the 
classical form of topological conjugacy, i.e.*? = ho&oh^ 1 for some homeomorphism 
h. We will show that they are equivalent in a very similar sense in the following 
section. 

3.2. CT-Equivalence. Let us find a reverse transformation from (|1.3[) to (|1.2p . It 
should transform f (^) back to r (£) which defines reverse timeshifts 

rjj = fj — rij, fj = max m . 

Analogously to (|3.ip , the reverse transformation T : C — > C is given as 
,„ ox ^ r 1 _ f (Vo)), (0) , * € [- min {%,r,} , 0] , 

l(yo(^ + %))j (o) , te [-^,-%-]. 

Again, we have 

(3.4) f of =$oT, 




Figure 3.1. CT-Transfo for a pair of bidirectionally coupled systems. 

and furthermore, T and T are reverse in the sense that 
(3.5) ToT = $^ and T of = % 

The equations f|3 . 2[) . (|3.4p and (|3.5[) constitute a notion of equivalence which can 
be applied not only to semidynamical systems which stem from systems of DDEs, 
but as well from delay coupled PDEs for instance. Therefore we define in a more 
abstract fashion. 
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Definition 2. Two semidynamical systems <f> : [0, oo) x X — >• X and : [0, oo) x 
Y — > Y are called CT- equivalent if there exist mappings TiX^y.TiF— >X 
and Tj > such that 

foT= fo* t = $ t of . 

We still use term CT-equivalence for the general formulation since, most proba- 
bly, the equivalent systems $ and ^ will in practice be connected by a transforma- 
tion as the componentwise timeshift (|3.1[) . 

3.3. Dynamical Invariants of CT-Equivalent Semidynamical Systems. In 

this section, we will infer some properties of CT-equivalent systems. Their state 
spaces hold a structure of corresponding strongly invariant sets. We prove, that sta- 
bility properties are preserved as well, presuming that the reverse CT-transformations 
T and T from Definition [2] are Lipschitz continuous. For the DDE systems 
and (|1.3|) this is the case if all fj are Lipschitz continuous. Throughout this section, 
we assume that X and Y are Banach spaces. 

Before stating the main theorem of this section, we give some definitions: 

Definition 3. A set A C X is called invariant if for any t > 0, 

$t (-4) = A. 

For any set A C X, we define its strongly invariant hull H§ (A) as the set 
i/$ (A) := {x S X | 3*i, i 2 > 0, & S A : $ tl (as) = $ ta (a;)} • 

A strongly invariant set A C A is a set that coincides with its strongly invariant 
hull, i.e. 

A = H$ ( A) . 

The class of strongly invariant sets is denoted by 

sis ($) = {ylCI | A = iJ$ (A)} . 

Theorem 4. Let $ : [0, oo) x I -> I and * : [0, oo) x Y~ ->• Y be CT-equivalent. 
Then, 

(i) for each (positively) invariant set A £ C, the set T [A] £ Y is (positively) 
invariant, 

(it) there is a one-to-one correspondence between (strongly) invariant sets. 

Proof. Ad (i): Let A be positively invariant, x £ A and y = T [x]. Then, ty t (y) = 
T [$< (a;)] G Hence T[A] is positively invariant. If A is invariant, then for 

each x £ A and each t £ [0, oo) there is an X-t £ A such that $t (x-i) = x and 
correspondingly for each y = T [x] £ T [A] there is y_ t = T [x-t] with 

9 t (y_ t ) =T[* t (a:_ t )] = T [a;] = y. 

Ad (ii): Let A G sis ($). We define a corresponding set A = Hy (T [A]) £ sis ($). 
Claim (i) is proven via 

o T) o o T) = id sis( $) and (Hy o T) o o T) = id sis (*), 
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where, by symmetry, it suffices to show only one equality. Let y E A, then there 
exist ti, t 2 > and y = T [x] E T [A] such that * t2 (y) = * tl (y) 6 A. Thus, 



f [* ta (y)]=T [* tl (y)] 

f [* tl (T[x])]=T [T[$ tl (x)]] 
$t 1+fJ (x) g ff* (A) = A. 



Hence, f [y] G A. That is, f 



that for a; G A, (x) £ T 



C A and ff$ 



(f[A\) 



C A. We have also shown 



and therefore x E ff$ (t A ^ . This yields 

A<ZH<s,(f A ) C A, i.e. A = ff$ A ) and (ff$ ofjo (if* o T) = id sis((&) . 
For invariant sets, the one-to-one correspondence is mediated directly via T and 



T. 



□ 



Definition 5. The maximal Lyapunov exponent (MLE) of a point x E X with 
respect to a semidynamical system $ : [0, oo) x I -) I on a Banach space X is 
defined as 



A (x) :— lim sup lim sup - In 



1, /|$*(x + 0-*t(aj) 



l£IV> 



G [0, oo] 



The MLE of a; G A C X with respect to the set A is defined as 



A (x, A) := lim sup lim sup min - In 



1 /|* t (x + 0-a| 



1 /dist($ t (x + 0,A) 



lei 



= lim sup lim sup — In , 

t->oo |£|\,0 * V 

The MLE of a set A G X is defined as A (A) = sup xeA A (x, A) . 

Theorem 6. Let the CT-trans formations T and T be Lips chitz- continuous and let 
X , Y be Banach spaces. Then, 

(i) corresponding (strongly) invariant sets of & and 'J possess the same maximal 
Lyapunov exponents (MLEs), 

(ii) for each positively invariant set A E C, the set T[A] E C has the same type 
of stability. 

Claim (i) is a direct consequence of the following Lemma. 

Lemma 7. Let X , Y be Banach spaces and let T and T be Lips chitz- continuous 
with constants Lt and Lf. Then, we have A(x) < A(T[x]) < A (x)) for all 
XEX. 

Proof. For each x, x G C, t > fj: 



$ t (x) -$i(jE + x)l 




(f o T [x]) - <i> t - n (f o T [x + x\) 



f o * t _^ (T [x]) - f o (T [a;] + £) 
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with |€| = |T[a; + x]-r[a;]|<Lr|xl 



< I i n ( L t ( T W) - *t-ij (T [s] + 01 



* V (1*1 /ir) 

1 ln ^^(T M )-^(T M+ 01^ + 1 ^ ^ 



Therefore, 



t \ \x\ J-t \ III 

and, thus, 

\{x) < A (T [a?]). 
The same reasoning for y = T [x] gives 

A(x) <A(T[x])< A ($,,(*)). 



□ 



Proof, (of Theorem E]) Lemma implies claim (i). To see that consider two corre- 
sponding sets A G sis ($) and A G sis (\&). Then it is impossible that A (A) > A (^A^j 
since for any x € A there exists y = T [x] E A such that A (x) < A (y). The same 
holds vice versa and therefore A (A) = A ( A ) . Similarly, one shows this for the case 



of invariant sets. 

Ad (ii): Let A G C be an (asymptotically) stable positively invariant set. We 
show that if A is (asymptotically) stable so is T [A). Let e > 0, y = T [xq] G T [A] 
and £ G C a small initial perturbation, i.e. |£ | < (5 with 5 = 6 (e) > to be 
specified later. Define the perturbed solution 

Vt = % (Vo + €o) = Vt + Ct. 
where y t = ^ t (y ) is the unperturbed solution. Consider 

2 t : = T [y t ] =f [*t(Vo+€o)l 

= *t(r[ y „ + €„]) 

= **(ror[aj ]+(r[ Vo + € ]-f [ Vo ] 
= + x, 7 ) , 



with 



X?7 



< Lrp ■ 5, 

where Lf is the Lipschitz constant of T. Since Xfj G A and A is (asymptically) 
stable, we can find a 5 = S (s) such that 

d(x t ,A) < 

and, in case of asymptotic stability, such that 

x t A. 



PARAMETER SPACE DIMENSION FOR NETWORKS WITH TRANSMISSION DELAYS 8 



This means, we can represent Xt as Xt — at + Xt with at G A, \xt\ < t~ an d, in 
case of asymptotic stability, \xt \ — > 0, for t — > oo. Note that t ^ at has not to be 
a solution. Define b t — T [a t ] E T[A]. Then, 



|» t+fl -6t| - ToT[y t ]-T[a t ] 
= \T[a t + Xt ]-T[at 
< L T \ Xt \<e 



and \y t+f j — 6*| —> 0, if A is asymptotically stable. This completes the proof. □ 



4. Reduction of Delay-Parameters 

The Theorems U and [5] show that CT-equivalence is indeed a very strong equiv- 
alence. For the most cases one is very well advised to study the systems which 
possess the most convenient distribution of delays within the concerned equiva- 
lence class. In section 2j] we will show, that it is always possible to find timeshifts 
rjj, j — 1, ...,N, such that the number of different delays in system reduces 
to the cycle space dimension C — L — N+loi the network, where L = is the 
number of links and N = #J\f is the number of cells. Effectively, this means that no 
more than C delay-parameters have to be taken into account during investigation. 
That this number cannot be reduced further is shown in section T4. 21 



4.1. Finding an Instantaneous Spanning Tree. 

Definition 8. A semicycle c— (£i,...,£k) is a closed path in the undirected graph 
which is obtained by dropping the orientation from all links. The delay sum of a 
semicycle c— (£i, ...,£k) with respect to a delay distribution is 

fc 

S r (c) :=^CTjT(^-), 

3=1 

where Oj E {±1} indicates whether the link £j points in the same direction as l\ 
(<jj = 1) or not (ctj = 0). The roundtrip of c is the modulus of its delay sum 

rt(c) :=|S(c)|. 

Lemma 9. Let r, f : £ — > [0, oo) be two delay distributions in a connected network. 
Then, the following statements are equivalent. 

(i) There exist rjj E K, j = 1, N , such that f {£) = t (£) — T]t(£) + VsU)- 

(ii) For all semicycles c of the network, 



S r (c) = S f (c) 
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Proof, (i) =>• (ii) Indeed, for any cycle c = (£i, ...,£k) we have 

k 

Sf(c) = 5> f &) 

fe 

= ( T (^) + %(/,)) 

k 

= S r (c) + a 3 (^%) _ ^ttf)) 

i=i 

= S T (C). 

(ii) (i) Select an arbitrary spanning tree S — (£\, £n-i)- Select £ = 
(£i,-,&v) and x = (Xi,-,Xiv) such that t (£j) = - and f(^) = 

Xt(^) — Xs(ij) f° r J = 1) •••) N — 1. (Note that both defining sets of equations are 
inhomogenous linear systems of type R JVx ( Ar ~ 1 ) and of rank N — 1. Therefore they 
possess solutions.) The shifts £ and % define distributions (with possibly negative 
values) 

h(£) = T (*)-&(/) 

f 2 (^) = f(^) -Xt(£) + Xs(£)- 

Both, fi and f 2 , are instantaneous along S, i.e., fi (^) = f 2 (•£) = for £ E S. Each 
link £ which is not in S (i.e., a fundamental link) corresponds to a semicycle c and 
is the only link in c which may hold a non-zero value fi (■£) or f 2 (£), respectively. 
By the part (i) => (ii) we have 

h (*) = S T (c) = S f (c) - f 2 (*) . 

Hence fi = f 2 . With i]j := £j — Xj this gives 

t(£) = T {£) - T] t(l) + T] s(Cj . 

□ 

Note that nj can always be chosen to be non-negative since a simultaneous shift 
rjj i-> ?7j + rj by some amount 77 e R does not change the resulting transformed 
distribution f . 

Let a network with delayed connections be given as a set of nodes TV = {1, TV}, 
a set of links £ C TV x TV, and a function r : C — > R>o which assigns to a link 
£ £ £ its delay t (^). For the link ^ = (k,j) which connects from node k to node j 
we define its source s (£) := k and target t(£) := j. The componentwise timeshift 
transformation [Eq. (3) of the main text] with timeshifts r]j, j = 1, ...,N, then acts 
on the delays as 

r {£) H- f (£) = r (f) + r7 s(£) - T] t{i) . 

Theorem 10. For every connected network with dynamics given by Eq. (1) of 
the main text, there exists a spanning tree S and timeshifts rjj such that in the 
transformed system (4) all links £ G S are instantaneous, 



f {£) = 0, for £eS, 
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Figure 4.1. Illustration of the fundamental section corresponding 
to the link £* = t\ +1 in the spanning tree Sk (solid links). Links 
which are not contained in Sk are indicated by dashed lines. 



and for each link I outside the spanning tree the delay f (£) equals to the roundtrip 
T(c{l)) [Eq. (5)j along the corresponding fundamental cycle c(£), 

f(£) = T (c (£)) , for £iS. 

Proof. The main idea of the algorithm is as follows. The timeshifts are changes 
interatively with r\ k denoting their values at the fc-th stage (k = 0, N — 1) such 
that, at stage k all links £ possess non-negative delays r k {£) > in the transformed 
system (4) with t]j — Vj- At each stage we create a spanning tree Sk containing 
k instantaneous links. Finally, we obtain timeshifts jj- — 1 and a spanning tree 

S = Sn-i whose links are all instantaneous in (4) with rjj = f]j~ l . Let us describe 
the algorithm in detail. 
Initial stage: 

Select a spanning tree So = {£%, ....,£ ( j v _ 1 } in the following way. Firstly, pick a 
link £\ with minimal delay, i.e., r (£j) = mmg e c { T (■£)}• Proceed adding the links 
£2, ■■■■,£ j with minimal delays as long as the set Sq — {£®, ...,1®} remains cycle- 
free. If at the j-th step the minimal link would create a cycle, do not add it but 
ignore it for the rest of this stage. We obtain the spanning tree Sq = Sq -1 . Let 
h E {0, ...,N — 1} be the number of instantaneous links, £ E So with t (£) = 0. 
Define S h := S , Vj °> = h -> N , and rh W = T W for ^ ^ £. Continue with 
stage h > 0. 

Stage k (k g {0, N — 2}): 

We are provided with timeshifts T)j , j = 1,...,N, such that the transformed 
system (4) with rjj = rf^ possesses a spanning tree Sk = {£*, ....,£%_ 1 } with exactly 
k instantaneous links. This spanning tree was constructed with respect to r fc {£) = 
t (£) + Vsfe) ~ Vttt) as ^0 was constructed in the initial stage with respect to r (£) . 
Choosing indices according to r fc < r k (£j +1 ), we have r k (ijj = 0, for j 

1 /:' and the link £* — £\ +1 has a minimal positive delay r k (£*) of all links in 

Sk- Take the fundamental section corresponding to £*. This means, we create the 
unique subsets of nodes Vk and Wk of N with s (£*) <G Vk and t (£*) S Wk, such 
that all other £ S Sk do not connect Vk and Wk (see Fig. 14. ip . 

Note that, by construction of Sk, the link £* = £ k +1 has minimal delay of all links 
between Vk and Wk- Since, if there were a smaller delay r k (£) < r k (£*) on another 
connecting link £, we would have picked £ instead of £* during the construction of 
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<Sfc. Next, we define new timeshifts r) k as 

~k = U k i+T k {t), forjGWfc, 
V 3 [Vj, ioi-jeV k . 
The transformed system (4) with rjj = fj k then has delays 



f k 



k (£) , for {s (£) , t (£)} C V k or {s {£) , t (£)} C W k , 

{£) = { r k {£) - r k {£*) , for s {£) G Vf. and t (£) G Wk , 



r k {£) + r k {£*) , for s {£) G W k and t {£) G F fc . 

We have f k (£ k ) = 0, for aU j = 1, k + 1, and f fc (^) > for all links £. Now, 
as in the initial stage, select a new spanning tree Sk+h beginning by adding the 
instantaneous links £ k , ...,£ k +1 . Then it has k + h > k instantaneous connections 
in (4) with r/j = J] k+h = fj k ■ If k + h = N — 1, the algorithm succeeded, elsewise, 
proceed with stage k + h. 

The statement that f (£) = T (c(£)) for £ ^ S 1 follows from Lemma [HI since 
f(T)>0. □ 

4.2. Genericity of the Dimension of the Delay Parameter Space. We call 
C the essential number of delays since it is the minimal number of different delays 
to which the number of delays in a network can be reduced generically. Here, 
generically means that the conditions which allow for further reduction of delays 
form a nullset in the parameter space of delays R> = {(r {£)) l££ '■ t (£) > 0} of 
the original system. The reducibility condition to m different delays is described 
by the L linear equations 

(4.1) ~ Vs W +t(1)=t(£), £e I,, 

with the restriction for f (£) to take one of m different values {6\, ...,6 m }. System 
(|4.1|) can be equivalently written in the vector form G 9 v = x, where v is the 
(N + m)-dimensional vector of unknowns v = (771, . . . , tjn, 9\ 1 . . . m ) and t is the 
i-dimensional vector of delays t{£). For any fixed assignment f (£) — 9 q (e), £ G 
Ij, with q : Ij {0,...,m}, and 6*o := 0, one obtains a different matrix G q G 
R L x ( JV + m ) . It can further be shown that the rank of the matrix G q is smaller than 
TV - 1 + rn. This implies for m < C = L - (N - 1) that N - 1 + m < L. Hence, 
the number of equations in (|4. 1[) is larger than the number of unknowns. Such 
equation cannot be solved generically, unless the given delays r (£) satisfy some 
special condition of positive codimension. 



